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Flavor changing neutral current (FCNC) interactions, 
which mediate transitions from one type of quark to 
another with the same electric charge, are suppressed 
in the standard model of particle physics (SM). FCNC 
processes are therefore sensitive indicators of physics 
beyond the SM (BSM). Presently there are only loose 
experimental bounds on FCNC decays of the t (top) 
quark [lj, the heaviest known quark. While the SM 
branching fraction for t — > Zq (q = u, c) is predicted to be 
0(1CP 14 ) 0, BSM models such as supersymmetry and 
quark compositeness allow branching fractions as high as 
0(1CT 4 ) @, 9. An observation of the top FCNC decay 
t —t Zq with present-size data samples would be a strong 
indication of BSM physics. 

In Tevatron Run I, a CDF search for the t — > Zq decay 
yielded the branching fraction upper limit B(t —> Zq) < 
33% (95% C.L.) 0. The current best 95% C.L. upper 
limit on B(t — > Zq), 13.7%, was set by the L3 experi- 
ment [5j from the non-observation of FCNC single t quark 
production. These analyses assumed a t quark mass of 
nit = 175 GeV/c 2 , close to the current world average of 
m t = 172.4 ± 1.2 GeV/c 2 0. In this Letter we also as- 
sume m t = 175 GeV/c 2 and evaluate the small effect of 
different nit on the result. 

We present the first Tevatron Run II search for the 
top FCNC decay t — ► Zq in tt pairs. The data sample 
analyzed for this Letter corresponds to an integrated lu- 
minosity of 1.9 fb _1 of pp collisions at y/s = 1.96 TeV, 
collected by the CDF II detector from March 2002 to 
May 2007. Our primary signal signature consists of the 
FCNC decay t — > Zq (q = u, c) together with the dom- 
inant SM decay t — > W~b (charge conjugate modes are 
implied). We search for the Z boson via its decays 
Z — > e + e~ and Z — ► and for the decay of the 

W boson into a quark-antiquark pair. Two leptons and 
at least four jets of hadrons, coming from the four quarks 
from secondary decays, can be observed in the detector. 
We also allow for events in which both t quarks undergo 
FCNC decays. We consider several SM background pro- 
cesses that also result in final states with a reconstructed 
Z boson and four or more jets. The dominant back- 
ground process is the production of Z bosons with as- 
sociated jets (-Z+jets). Smaller contributions come from 
SM tt production and the production of pairs of gauge 
bosons (dibosons), WZ and ZZ. The contributions from 
WW diboson production and from W bosons produced 
in association with jets are negligible. 

The components of the CDF II detector relevant to this 
analysis are briefly described here: a more complete de- 
scription can be found elsewhere Ml. The transverse mo- 
menta px and pseudorapidities r\ |8j of charged particles 
are measured by a silicon strip detector Q and a 96-layer 
drift chamber (COT) [l(| inside a 1.4 T solenoidal mag- 
netic field. The silicon detector and the COT provide 
good combined reconstruction efficiency for \rj\ < 1.2. 
The precise position resolution of the silicon detector is 
also crucial for identifying displaced secondary vertices 
from from long-lived B hadrons (6-tagging) . Electromag- 
netic (EM) and hadronic calorimeters measure energies 



of charged and neutral particles in the central (\n\ < 1.1) 
and end-plug (1.1 < \n\ < 3.6) regions. Drift cham- 
bers and scintillation counters located on the rear of the 
calorimeters and behind an additional steel absorber de- 
tect muons with \r]\ < 1.0. A three level trigger selects 
events that contain electrons (muons) with Et > 18 GeV 
(p T > 18 GeV/c). 

We use the pythia v6.216 Monte Carlo (MC) genera- 
tor [ll| to simulate the FCNC signal in tt events and 
all sources of SM background, except Z+jets produc- 
tion, whose kinematic distributions we simulate with the 
ALPGEN MC generator [12], v2.10' interfaced to PYTHIA 
v6.325. Acceptance, efficiency, and kinematic distribu- 
tions of signal and background are determined from the 
above MC simulations. We add a t — * Zq decay chan- 
nel to PYTHIA by forcing the Z helicity to be consistent 
with the expectation from an SM-like Higgs mechanism 
of 65% longitudinally polarized and 35% left-handed Z 
bosons. The main signal MC sample contains events that 
decay according to our primary signature, tt — > Zq Wb 
with q — c. If the t — > Zc decay is replaced by t — > Zu, 
the probability to 6-tag an FCNC signal event is re- 
duced from 50% to 45%, where we apply the SecVtx 
algorithm Q for the 6-tagging. We have also gener- 
ated an MC event sample in which both t quarks decay 
via the FCNC mode and the Z is allowed to decay into 
the e + e _ , and qq decay modes. The amount of 

tt — > Zc Zc relative to tt — > Wb Zc is adjusted accord- 
ing to the branching fraction B(t — > Zq) during the limit 
calculation. 

We apply corrections for trigger efficiencies and for dif- 
ferent efficiencies in data and simulation for lepton iden- 
tification and reconstruction, 6-tagging, and misidentifi- 
cation of secondary vertices (mistags). 

The analysis proceeds in three steps. Our base selec- 
tion requires a Z boson and four or more jets. An op- 
timized event selection separates the FCNC signal from 
the SM backgrounds based on kinematic properties of the 
event. All events that pass the optimized selection are 
further divided into two signal regions based on whether 
any of the four leading jets has an identified fe-jct ("6- 
tagged region") or not ( "non-6-tagged region"); those 
events that pass the base selection but not the optimized 
selection are placed in a "control" region. While only 
12% of FCNC signal acceptance falls in the control re- 
gion, two-thirds of the SM backgrounds are located there. 
We finally fit the data to one-dimensional probability 
density histograms ( "templates" ) of one of the kinematic 
distributions (mass y 2 , see below) to derive a limit on 
B(t — > Zq). This is illustrated in Fig. Q] To reduce sys- 
tematic uncertainties, the measurement is normalized to 
the measured event yield in SM tt production. 

Our base selection chooses events with two oppositely- 
charged leptons of the same flavor (e or fi) and four or 
more jets. One lepton is required to be a central electron 
or muon, the other can be a forward electron or a track, as 
described below. The leptons must be compatible with 
originating from a Z boson in the mass window from 
76 GeV/c 2 to 106 GeV/c 2 (> 3cr in units of CDF's Z 
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mass resolution). 

All lcptons used in this analysis arc required to be well 
isolated in a cone of AR = ^(A?y) 2 + (A0) 2 < 0.4 @| 
and to have transverse energies Et (momenta Pt) greater 
than 20 GeV (20 GeV/c). Electrons are identified by re- 
quiring an energy cluster in the EM calorimeter with a 
single track pointing to it. Central electrons are required 
to have a high quality COT-based track, calorimeter clus- 
ter Et consistent with the track pt, a high fraction of 
the total energy deposition in the EM calorimeter, and 
a lateral shower profile consistent with electron show- 
ers. Forward electrons are reconstructed in the end-plug 
calorimeter and have similar constraints except that the 
tracks are reconstructed only in the silicon detector and 
the calorimeter cluster Et and track pt are not com- 
pared. Muons are identified by matching tracks recon- 
structed in the COT to track segments reconstructed 
in the muon chambers and by requiring energy deposi- 
tions in the calorimeters consistent with minimum ion- 
izing particles. The muons in this analysis are required 
to be in the central region (j^l < 1.0). We double the 
acceptance for leptonic Z decays by allowing one of the 
two lepton candidates to satisfy weaker selection criteria, 
requiring only an isolated track which passes COT and 
silicon detector track quality cuts. For tracks used as 
electrons, if an EM calorimeter tower is associated to the 
track and the energy of the EM tower is greater than the 
track momentum, the EM tower energy is used instead 
of the track momentum. 

Jets are identified by energy deposited in the calorime- 
ters within a cone of AR < 0.4. To improve the parton 
energ y e stimate, jets are corrected for instrumental ef- 
fects [13[. We select events with at least four jets with 
corrected E T > 15 GeV and | ^7 1 < 2.4. 

We separate the FCNC signal from the MC back- 
ground with the help of further selection criteria in ad- 
dition to the above base selection that will form the op- 
timized selection, a mass x 2 , the transverse mass of the 
system, and the Et of the four leading (highest Et) jets. 

The decay ti — > WbZq — > qq'bttq" contains no high- 
energy neutrinos; therefore we can fully reconstruct the 
event kinematics. The four jets in signal events result 
from the b quark and the decay products of the W in 
the t — > Wb decay, and the c or u quark from the t — > Zq 
decay. We form all permutations of the four leading jets 
in the events to compare the reconstructed masses (m rec ) 
of the W, top quark decaying to Wb, and top quark 
decaying to Zq. We define a mass x 2 as 

2 _ ( m W,rcc — "Iff \ / m t ^Wb,iec ~ ™t \ ^ 

\ &W J \ &t-,Wb ) 

+ / mt^zg.icc - m t \ 2 

V Pt^Zq J 

and select the permutation with the smallest \ 2 . We 
scale the measured four-momenta of the W and Z boson 
daughter particles such that the boson masses are fixed 
to the world average values [3] and use the scaled four- 



momenta to calculate the two top quark masses. The 
widths used are given by the standard deviations of the 
reconstructed masses measured in the MC simulation of 
FCNC events. Using the correct pairing of jets to par- 
tons, we extract aw = 15 GeV/c 2 , a t ->wb = 24 GeV/c 2 , 
and <j t ->zq = 21 GeV/c 2 . We expect FCNC signal events 
to populate the low x 2 region and background events to 
result in higher x 2 , see Fig. [TJ We have verified that the 
components of Eq. ([1]) describe the data well in events 
with a Z boson and three jets. 

Since the FCNC signal events originate from ti de- 
cays, they contain more central Z bosons and jets than 
background events. To exploit this, we use the trans- 
verse mass of the Z and the four leading jets, defined 

as ttit = \J (Y^, Et) 2 — (J2pt) 2 , as a selection criterion. 
We also apply a tiered cut on the Et of the four lead- 
ing jets, as FCNC signal events contain jets with higher 
transverse momenta than SM background events. 

We optimized these additional selection criteria for the 
best expected limit on Bit — > Zq) in the absence of a sig- 
nal, using the MC simulation and a signal-depleted con- 
trol region in the data (\fx* > 3). The optimization 
was performed in the blind phase of a counting experi- 
ment analysis using the first 1.1 fb _1 of integrated lumi- 
nosity [HI]. We leave the optimized selection unchanged 
for the full 1.9 fb _1 result; it requires transverse mass 
m T > 200 GeV/c 2 , leading jet E T > 40 GeV, second jet 
E T > 30 GeV, third jet E T > 20 GeV, and fourth jet 
E T > 15 GeV. After optimization, 88% of the FCNC sig- 
nal events from the base selection fall into the two signal 
regions, compared to 33% of the background events. The 
inclusive signal acceptances for the decay ti — > WbZc 
(ti -»■ ZcZc) are 0.43% (0.58%) for the 6-tagged selec- 
tion, 0.34% (0.86%) for the non-6-tagged selection, and 
0.10% (0.16%) for the control region. 

To determine the FCNC branching fraction, we take 
into account single or double FCNC decays of ti pairs 
and normalize to the event yield of a selection for the SM 
decay ti — > WbWb — > Hvbqq'b ("lepton+jets") requiring 
at least two jets to be secondary vertex 6-tagged [16[ . In 
1.9 fb _1 we observe 277 ti candidate events, consistent 
with a production cross section of 8.8 ± 0.7 (stat.) pb 
assuming B(t — > Wb) = 100%. If t — > Zq decays were 
present, these additional ti decays are less likely to be re- 
constructed in the lepton+jets mode, resulting in a mea- 
sured ti production cross section smaller than the actual 
cross section. We correct for this effect by modifying 
the measured cross section based on the limit we set on 
B(t -► Zq). 

We extract a limit on the branching fraction Bit — > Zq) 
from a fit to the mass x 2 distribution using templates 
constructed from the MC simulated mass x 2 distributions 
of the FCNC signal and the SM backgrounds (Z+jets, 
SM ti, and dibosons). The normalization of the dominant 
Z+jets background is the most difficult to estimate from 
data and MC simulations; therefore it is extracted from 
the fit. The SM ti background is normalized to the ob- 
served event yield in the lepton+jets decay mode; back- 
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TABLE I: Results of the fit to data. From the ratio of the 
number of Z+]ets events in the signal and the control regions 
and the tagging fraction /tag, we obtain Ztaggcd = 13.5 events 
and ^non-tagged = 53.9 events. 



Fit Parameter 


Value 




Branching Fraction, B(t — > Zq) (%) 


-1.49 


± 


1.52 


Z+Jets Events in Control Region, Z coniIO \ 


129.0 


± 


11.1 


Ratio Signal/Control Region, 7£ s ig 


0.52 


± 


0.07 


Tagging Fraction, / tag 


0.20 


± 


0.06 


Jet Energy Scale Shift, erjES 


-0.74 


± 


0.43 



grounds from diboson production are normalized to their 
theoretical cross sections. Both contributions arc fixed in 
the fit. The expected background from SM ti production 
and diboson production is 2.2 ± 0.2 (3.2 ± 0.2) events 
for the 6-tagged (non-6-tagged) selection. The 6-tagged, 
the non-6-tagged, and the control region are fit simulta- 
neously. We include systematic uncertainties due to the 
shapes of the signal and background templates by allow- 
ing the templates to change shape via a histogram inter- 
polation technique (horizontal template morphing) [TtJ • 
The uncertainty due to the experimental jet energy scale 
(JES) is more than three times larger than any other un- 
certainty, and the data are most consistent with a shift 
of the JES; therefore we only include JES induced shape 
uncertainties in the fitting procedure. 

We use the signal-depleted control region to constrain 
the background shape uncertainties without losing sen- 
sitivity to a small FCNC signal. Additionally we use 
the number of Z+jets events observed in the control re- 
gion, Z contro i, to place a loose constraint on the number 
of Z+jets events in the two signal regions, Signal- We 
constrain the ratio 1Z S i g = Signal /^control to the value 
estimated by the MC simulation, 7\L S i g = 0.51 ± 0.10 at 
the nominal JES. The uncertainty is conservatively esti- 
mated by varying the energy scales in the ALPGEN MC 
generator. We adjust 7?. s i g as a function of the JES shift 
ctjes, keeping the relative uncertainty of 20%. The ab- 
solute number of Z+jets background events remains un- 
constrained. To reflect the constraint on 7£. s j g in the tem- 
plate fit, we parameterize the number of .Z+jets events 
passing the 6-tagged and non-6-tagged signal selections, 

^tagged and ^non-tagged , ^tagged ./tag ' 7^-sig ' ^control 
and ^non-tagged = (1 — /tag) 1 7?-sig 1 ■Zcontrol, where /tag 

is the fraction of Z+jets events passing the 6-tagged sig- 
nal selection. In summary, the unconstrained parame- 
ters of the template fit include the branching fraction 
B(t — > Zq), Z CO ntroi, the tagging fraction / tag , and the 
shift in JES <t.jes • We apply a Gaussian constraint on the 
ratio lZ s i g . From the template fit to the data we measure 
a branching fraction B(t — > Zq) = —1.49%. The fit result 
is summarized in Table U and in Fig.[TJ 

We employ a Feldman- Cousins (FC) limit calculation 
framework that includes the handling of systematic un- 
certainties [HI- The FC construction is based on ap- 
plying the above template fits to simulated experiments 
that are generated taking into account all known sources 
of systematic rate and shape uncertainty and their cor- 
relations (systematic uncertainties are discussed in the 
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FIG. 1: Mass \ 2 distribution for fr-tagged and non-&-tagged 
signal regions and the control region. The data points as a 
function of \l~x* are compared to the SM background predic- 
tion and the expected FCNC yield at the observed 95% C.L. 
upper limit on the branching fraction B(t — > Zq) < 3.7%. The 
data are consistent with the background prediction. 

next paragraphs). We expect a limit of (5.0 ± 2.2)% in 
the absence of signal. From the FC construction we find 
B(t -> Zq) < 3.7% at 95% C.L. 

We have studied systematic uncertainties due to shape 
uncertainties of the fitted templates, rate uncertainties 
(signal acceptance, background rate), and due to the 
normalization to the lepton+jets event yield. We treat 
the dominant source of template shape uncertainties, 
the experimental JES, MS M. free fit parameter. To ac- 
count for further shape uncertainties, we measure the 
bias on the fitted branching fraction B(t — > Zq) in simu- 
lated experiments generated assuming the second largest 
source of shape uncertainties, the MC generator used for 
the iv+jets background. We simulate Z+jets events for 
which the vertex energy scale is varied by factors of two 
and assign the measured relative branching fraction bias 
at B(t — ► Zq) = 3.7% as a systematic uncertainty (5.6%). 

The rate uncertainties, summarized in Table [TT] are 
divided into anti-correlatcd uncertainties, those which 
cause migration of events between 6-tagged and non- 
6-tagged selections, and correlated uncertainties, those 
which simultaneously increase or decrease both selec- 
tions. The signal acceptance systematic uncertainties are 
evaluated for the ratio of the FCNC signal acceptance to 
the acceptance for the event selection used in the lep- 
ton+jets normalization mode. Background rate uncer- 
tainties affect only the smaller SM ti and diboson back- 
grounds. The rate of the dominant Z+jets background 
in the control region is a free parameter in the template 
fit; therefore we do not assign systematic uncertainties 
for this. 

Normalization to the lcpton+jets event yields in the 
ti production cross section analysis removes nearly all 
uncertainties depending directly on luminosity. Many 
other systematic uncertainties also partially cancel. Wc 
absorb the statistical and systematic uncertainties of the 
ti production cross section measurement as part of the 
systematic uncertainty (7.8%). 

In conclusion, we have searched for the top quark fla- 
vor changing neutral current decay t — > Zq in events with 
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TABLE II: Systematic uncertainties of the FCNC signal ac- 
ceptance relative to the acceptance of the lepton+jets normal- 
ization mode and the background rate for the 6-tagged and 
non-6-tagged selections. 



Systematic 


Acceptance 


Background 


Uncertainty 


b 


Non-b 


b 


Non-b 


ID, Trig, and PDF 


1.0 


1.0 


1.5 


1.5 


Initial/Final State Rad. 


4.8 


5.5 






Z Helicity 


3.4 


3.6 






Total Correlated 


6.2 


6.1 


6.2 


6.2 


^-Tagging 


5.6 


16.1 


3.2 


2.5 


B{t -> Zc) vs. B(t -> Zu) 


4.5 


4.5 






Total Anti-Correlated 


7.2 


16.7 


3.2 


2.5 



a Z boson and four or more jets using CDF Run II 
data corresponding to 1.9 fb _1 of integrated luminosity. 
The data are consistent with the SM background predic- 
tion, and we set an upper limit on the branching fraction 
B(t — > Zq) of 3.7% at 95% C.L. at a top quark mass of 
175 GeV/c 2 . Assuming a top quark mass of 170 GeV/c 2 , 
the 95% C.L. upper limit is 4.1%. Compared with the 
previous world's best limit of 13.7% by the L3 experi- 



ment [5j at LEP2, and the previous Tevatron upper limit 
of 33%, as reported by CDF Run I our reported limit 
of 3.7% represents a substantial improvement. 
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